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THE MIXTURE PROBLEM IN COMPUTER MAPPING OP TERRAIN: 

IMPROVED TECHNIQUES FOR ESTABLISHING SPECTRAL SIGNATURES, 

ATMOSPHERIC PATH RADIANCE, AND TRANSMITTANCE 
By Harry W. Smedes, U.S. Geological Survey, Denver, Colorado, 

Roland L. Hulstrom, Martin Marietta Aerospace Corporation, Denver, Colorado, and 
K. Jon Ranson, Colorado State University, Fort Collins, Colorado 

ABSTRACT 

Among the results of our LANDSAT and Skylab research programs on the effects of the 
atmosphere on computer mapping of terrain, we emphasize the following: (l) The concept of 

a ground "truth" map needs to be drastically revised; (2) the concept of training areas and 
test areas is not as simple as generally thought because of the problan of pixels that re- 
present a mixture of terrain classes; (3) this mixture problem needs to be more widely re- 
cognized and dealt with by techniques of calculating spectral signatures of mixed classes, 
such as those we used, or by other methods; (U) atmospheric effects should be considered in 
computer mapping of terrain and in monitoring changes; and (5) terrain features may be used 
as calibration panels on tte ground, from which atmospheric conditions can be determined 
and monitored. 

The natural terrain comprises a mixture of diverse classes, including combinations of 
such things as tree, grass, and other vegetative- cover species, different types of soil and 
rock, and water bodies of different size, depth, and clarity. Because the terrain features 
generally are small compared to the ground resolution element (= pixel) from multispectral 
scanner data, it is unusual for very many contiguous pixels to consist of a single terrain 
class. Commonly, each pixel is a mixture of two or more classes. 

It has been shown that the spectral signature (radiance) of a mixture of terrain 
classes is not representative of any of the component classes composing the mixture. Auto- 
matic recognition processors may therefore misclassify mixture- containing pixels, which 
often leads to an underestimation of the amount of component classes present in the area 
scanned. 

Even a single class may have a range in spectral signature depending on such things as 
size of pixel (sample size), solar aspect, orientation, and atmospheric conditions. 

Because of this problem of mixtures in pixels, we at first attempted to compile ground 
control maps whose classes were designated by different increments of the more 
ring mixtures and to train the computer to recognize th«n by selecting specific TRAINING 
AREAS of known proportions of classes. Not only was this a monumental task, loaded with 
subjective judgements and difficult to calibrate, but the TESTING of accuracy of the result- 
ing computer-derived map proved to be equally monumental. After much effort we concluded 
that— if not lmpossible--it was certainly not at all feasible to prepare a truly accurate 
ground control map. In fact, while checking for errors in the computer map, the ground 
control or so-called ground "truth" map continually had to be upgraded. A serious problem 
that needs to be faced is: '"what constitutes the ground truth?' Our first major conclusion 

was that, if properly trained on end-member and mixed classes, the computer made a more 
accurate map of terrain cover than we were able to compile as control data. 

However, it is very difficult to measure precise proportions of mixtures and to locate 
those areas in terms of LANDSAT MSS, Skylab S-192, or other scanner pixels, especially in 
areas of mountainous terrain. This difficulty and the resulting uncertainty and imprecision 
led us to experiment with techniques to calculate mixes from data for end members and to 
use those calculated signatures in lieu of training areas for the mix classes. Previous re- 
search in this problem has largely been confined to classes of crops in flat agricultural 
fields. Our data and techniques will apply equally well to conditions ranging from crop- 
lands to wilderness. 


Inasmuch as the radiance from a pixel is integrated over the entire area of the pixel 
Instantaneous field of view of the scanner), mixture- class signatures can be deterained' 
^ the known mean vectors and covariance matrices of the component classes and thepropor- 

Signatures for mixed classes w^ 

alculated using spectral radiance from on-site measurements. In addition, several experi- 
ments were conducted with simulated LANDSAT MSS-type data to show the exnected Inmmw^pnt 

classification using simSated miiure'liL* s1f^Ll2e.?™o 
determining component class signatures when there were insufficient 
data points for conventional signature extraction. 


for terrain of south-central Colorado 

toJefand mixture-class spectral signs- 


these techniques highly successful in terms of shorter preparation time 
and higher accuracy, but the calculated signatures also afford insights into what the appro- 

discrLlJISl^r^* mixes and what the optimum wavelength bands are for the most accurate 
discrimination among specific classes. 

Atmospheric effects must be known in order to combine spectral signatures derived 
d^tly from LANDSAT or Skylab scanner tape data with those measured on-site and from air- 
craft. A few large homogeneous sites such as a large body of deep clear water, a lante 
expanse of bare rock, dense forest, large nonvegetated dry lake bed or desert sand, and 
snw may serve as known natural calibration panels on the ground. These would be visible 

atmospheric path radiance and transmittance can be derived 
^ flight). This calculated path radiance and transmittance 

sense, this technique can serve as a useful means of monitoring atmospheric 

Sharps ^Ircrax't, as atmospheric path radiance and transmittance are 
measures of air quality. 


INTROIWCTION 

data from satellites are used as input to computers for automati- 
ground cover. Some major problems faced in this remote- 
sensing task include 1) the effect of mixtures of classes and, primarily because of mixtures 

sl^nJ^^rir constitutes ground "truth"; and 2) effects of the atmosphere on spectral* 

paper presents the fundamental principles of these problems and some of 
the results of our studies of them for a test site in Colorado, using LANDSAT-1 data. 

♦ — (figure 1) comprises about 2280 sq. km (880 sq. mi) of generally wildland 

Colorado. It includes such landmarks as Pikes Peak, the Cripple 
S hlchw^Cfn and Canon City. Altitudes range from 1525 to U300 m. The te?Sin 

iide^SiL ^ ? includes a diversity of rock types, soil, and vegetative cover, over a 
wide range of angle and aspect of slope. * 

funded'bv*‘SI^*’‘'!!n^KvT'’J®^ an integral part of LANDSAT and Skylab EREP projects 

by NASA, and by ^-house research supported by the Martin fterietta Aerospace Corpo- 

mn Si; m^df by Colorado State University and a computer-derived 
map was made by ^e Environmental Research Institute of Michigan (ERIM) as their part of a 

c matures problem. This research effort was conceived, initiated and 

^rdinat^ by Smedes. Hulstrom measured atmospheric properties and studied their effects, 
mode on-site measurements of spectral signatures, and used these data in a con9>uter prograii 

mixtures, tenson studied the effects of mixtures using simulated 
LANDSAT data as part of a dissertation for an advanced degree at Colorado State University. 

THE MIXTORE problem 

The natural terrain is composed of mixtures of classes. This is true regardless of the 
size of the ground resolution element ( ■ pixel) of the sensor system, and holds for satellite 


data with Resolution on the order of 100 sq. meters, to that of the microscope with reso- 
‘ lution of the order of a sq, micron or less. tradition, each discipline of the natural 
sciences has accommodated this mixtures problem by means of graded orders or hierarchies 
of classifications. We thus havt'. In order of decreasing size of resolution element, such 
classes as galaxies; stars and planets; continents and oceans; mountains, ridges, valleys, 
ptlains and deserts; woodlands, meadows, bare rock, and soil; spruce, aspen, bunchgrass, 
granite, and quarts veins; bark, needles, leaves stems, quarts, feldspar, and mica grains; 
petioles, chloroplasts, mitochondria; perthite, alblte with fluid inclusions , and magnetite 
inclusions in biotite. Because we are dealing with current satellite scanners, the reso- 
lution element or pixel size we are concerned with is 6^-8o m square. Tills resvdts in high 
likelihood that any given pixel will contain two or more of such classes as wtxHllands, 
meadow, bare rock, and soil. If a recognition processor (computer) is trained to recognise 
only the homogeneous component classes (trees, grass, rock, etc.) the overall consequences 
of a mixture of these classes occurring within a pixel is a misclassifled or unclasslf leti 
pixel I which tends to make estimates of the area covered by a terrain class lower than 
actually present (ref. 1). Errors in classifications due to the mixture problem alone will 
be 25 to 30 percent (ref. 2 and 3). 

There are two basic kinds of mixtures. One is due to the boundary or e<lge effect, 
where a pixel encompasses the interface of two or more homogeneous end-mianber component 
classes, such as water and dense coniferous forest. Tlie resulting spectral signature would 
not be representative of a pixel that containt'd only one of those classes. Ttte other mixture 
is representtKl by pixels which contain homogeneous mixtures of such things ns forest canopy 
and grass, rock, or soil understorj'. In the first case, only the pixels along the inter- 
face will bo misclassified. In the socond case, largo clusters of pixels throughout the 
forest will be misclassifled if otvly tl»e component classes were ustvi to train the comptiter 
(ref. U, S, 6). 

variations in the amount of vegetative cover will result in corresponding mixtures of 
the vegetation and underlying material. Consequently, if the ci>mputer is traineii to recog- 
nise only one vegetation density, then misclassifications would occur for densities above 
or below that of the training class, within the constraint of some rcsi>onsc threshold. These 
misclassifications occur because when a mixt.uro of terrain classes is contained within the 
instantaneous field of view (rf*0V) of n scanner, the spectral respoitse obtained is unlike 
that of any of the ^'oinpcnent class signatures (ref. 8). To illustrate this concept, three 
terrain classes were Identificvi on grajanaps of August 20, l')72 (frame no. 1028-1713*>) LANB- 
SAT-1 MSS data that appivinxl r«'presentntive of grassland (A), dense forest assumed to be 
100^ cover (B) and a class representing an assumed SO-SO mixture (C) of classes A and B 
and located on tiie interface between them. Spectral signatures were extracted with 21 data 
points sampltxl for Class A, 72 for class B and 2U for class C. The mean spectral response 
curves for the three classes arc shown in figure 2. Note that the spectral response for 
the mixture class C is uncharacteristic of either of its component classes (A and B) but 
that the response curves falls between the two iximponents' curves. Researchers investigating 
this phenomenon have found that the relationship of a mixture and the signatures of the 
component classes may be appn^ximatod as a Ainctlon of the area of the pixel occupied by 
each component terrain class and their respective spectral sigiiatures. 

Because of this problem of mlxi.ures in pixels, we at first attempted to compile ground 
control maps whose classes were designated by different incriwents of the more widely occur- 
ring mixtures and to train the computer to recognise them by selecting specific training areas 
of known proportions of classes. Tt’-ls methevi, described in detail in ref. 6 and 7, is ex- 
tremely time-consuming, loaded with sub.lectlve Judgtsnents, and very difficult to c^lbratc. 
Testing the accuracy of the resulting computer-deriveti map is even more tlme-i'onsvuning 

Ckir conclusion, supported by other work (such as reported in ref. 6, 7, and d), was 
that if the computer is properly tmineti on end-member components and on commotily-occurrity^ 
mixtures, it (x>uld moke a more accurate map of the terrain ettver than we were able to cewpile 
as control data. TlJin oasts serlovis doubt upon atul calls for reconsideration of what consti- 
tutes the ground 'Hrutli''. 

However, it la veiy difficult to measure precUa proportions of mixtures and to looata 


thoae areas in terms of LANDSAT MSS, Skylab S-192, or other scanner pixels, especially in 
areas of mountainous terrain. This difficulty and the resulting uncertainty and imprecision' 
led us to experiment with techniques to calculate mixes from data for end members and to 
use those calculated signatures in lieu of training areas for the mix classes. Previous 
research in this problem has largely been confined to classes of crops in flat agricultural 
fields. Our data and techniques will apply equally well to conditions ranging from crop- 
lands to wilderness. 

Inasmuch as the radiance from a pixel is integrated over the entire area of the pixel 
(the instanteineous field of view of the scanner), mixture-class signatures can be deterained 
from the known mean vectors and covariance matrices of the cooqponent classes emd the propor- 
tional areas occupied by each component in the pixel. 

TWO approaches were taken. One was to measure the reflectivity of pure end-member 
components on site and then calculate signatures of various mixtures. The other was to 
use the satellite data itself to extract the radiance signatures of known end-member targets 
and then calculate radiance signatures for various mixtures. Simulated LANDSAT MSS- type 
data were used. They show the expected improvement in accuracy of autcsnatic classification 
using simulated mixture-class signatures. In addition, methods have been developed for de- 
termining component class signatures from mixtures when there were insufficient data points 
for conventional signature extraction of the end-member components. 

CALCULATED SPECTRAL SIGNATURES FOR MIXTURES 


Basic Concept 


The basic problem of correctly identifying ground resolution elements that contain mix- 
tures of terrain classes with simulated mixture signatures is obtaining representative sig- 
natures for the component classes. A component class can be defined as a homogeneous or 
non-horoogeneous group of materials that make up a discrete mapping class or terrain type. 

A discrete mapping class may or may not be found in combination with other discrete classes 
as mixtures. For example, in an area where there is a sparse covering of grass with bare 
soil showing through, the grass and soil may be considered as two discrete mapping classes. 
With the O.U ha resolution of LANDSAT, the scene would be viewed as a mixture of grass and 
bare soil so these two classes would most likely be combined into one component class for 
analysis. If another class were identified as dense forest, then euiy time enough trees 
existed on the grass-soil unit to affect the response of a pixel, then a mixture of these 
two classes (forest and the grass-soil association) would exist. Discrete mapping class 
determinations can usually be considered a function of the natural associations, the objec- 
tives of the user, spectral radiance, and the limitations of the MSS data. 

When selecting component classes for mapping mixtures, care must be taken to avoid 
situations where the data for a class form a multimodal distribution. Multimodal distri- 
butions occur where the response data for a class are affected by variables such as slope 
and aspect, vegetation vigor, underlying soil spectra, and sensor scan angle effects. The 
standard method for dealing with this problem is to divide the multimodal class distribution 
into subclasses as a function of slope, etc., classify the data and then combine the results 
for each subclass (ref. 10, ll) for final display. 

Calculations From On-site Measurements 


Mixtures.- For on-site and near-surface measurements (negligible atmospheric 
and path radiance) the radiance, N^ of a pure target end-member can be given as 


N 


t 


HP 

TT 


transmission 

( 1 ) 


Where H is the total solar irradiance and P is the target reflectivity. Tlie radiance from 
a mixture of targets, Nm, can be given as 


N 


m 






U 


( 2 ) 


whertP . P , etc. are the reflectivity signatures of each component target, and 

etr aVi tfe fractional proportions of the pixel that t^ach component target 

«^4tton !, ■«<». that Son; of the component tarfots have atrong apocular rnn«ttyltl«. 

Note that even for the nteir-surfaee conditions, there are atmoapherio effecta to be 
considered (H, in equations 1 and 0, above). 

In onier to siimdate mixed target radiance signatures, a knowliHige of the total in- 
eominc solar irradiance (H), the pure target reflectivity signatures, and the fractional 
art>a of each target need to be known. In this study the spectral distribution and magni- 
tude of the Ini-oallnB aolar irradiance ia r.eaaurwl for v«rloua '‘''Kj;” ^ S hSuroSter 

S™rttih::i" arSS r:"pre.7S5 «“:ttU i» si^ply varied to ,i..a«te varioua oc«bihatlo„. 
of mixtures. 

The first mixture problem addresseti was that of Pikes Peak Granite and f®'** 

Forest This representl a very r.lmple and ccmnKMi mixture. Hie reflectivity 

Peak Granite and^he total Incident solar irradiance were measunnl on site (as discussed in 

tht^following sections) while the refleet ivity signature for Coniferous 

from ref Ict Figure i shows all three of the parameters, -nie resultant simulations of 

eomoonent and mixture signatur.'s for the LAUDilAT bands and for the continuous 

region from 500 to 1100 im are shown in figures U and 5. It is ® 

(as shown In^fig. U ami 5) for a ^0-50 mixture of rock and trees, the resultant mixture 
ridlance signature no longer contains the characteristic chlor>’phyll absorption 

^«^tlhrr«5Lco .^nature would neat Hkoly b. iht.-rpretod as a grayish. lov-alb«lo 

rock. 

<5noet.ral contrast - In addition to con\puter oimiaations of pure and mixed target radi- 
anoo gr.rt -^;^c-tra? Inhoront oohtnhst of various tarots and th.lr ^ok^round can 
be colcSlate<i. Tlu* inherent contrast. C, of a target versus its Ivickground is 

^ ( 3 ) 


N, 


Where N is the rnainnce signaun-e of the target and Is the r«dl«.H-e signature of U 
backgroi^mt. The target and/or background c«,» be any p?ire target or mixture ot targets. 

Smctral contrast was calculated In onier to determine if there existed any particular 
spectral region or band that gav, a better contrast U-.an other spectral 

The goal o/these spectral contrast calcnlat Ions Is not to exclusively attempt to dtttrmlnr 

if "ontlmum" I'ands existed for the dlscrlmlnat ion of rtn'k iyp« s, vegetation t.vpes, etc. 

tother, the givil is to determine "optlmiun" bands for maximum spectral contrast 

selected target(s) and Its background; because, for the natural terrain existing in the 

Colrrlido t«T.u’, tblk U in Hu., will, tu.- uctuu. t-rr.!.. coudllUulk 

two rock types rarely exist slde-by-slde with good expi^sure; instead, a hiv.n ro^ t>pt is 

«po,ed duil .urroun.!;.! by u l«..kpn..,nd n«d.. up of u of 

thrrock types do exist side-hy-side they usually do not till tiie entire ut oi 

the sensor. This is dne to the fact that they occur as small outcrops and/ or they exl-t 

layers exposed on extremely st»*ep slopes or clitts. 

An example of the computer sp«-ctral contrast calcidatlons is shown in figure 6. The 
pure target is Pikes Peak Gnmlte. and the Ivickground consists of mixtures of «'**“' ^** *"‘* 
Lntferous fort'st It is obvious that there an* optimum spectral regions; and, if tht 

crprl;es a mixture, the contrast Is reduced tremendously. The same case only 
aptxlled to the^broader bfinds of UNDSAT. Is shown In figure 7. As can be optimum 

^UriUever. the .werall contrast is reduced because of the br.v.der l>and.. 

Ttu' eonputer calculations of spectral contrast are not limited to the UNDlUT 
re-a^irtiorSlt ca^nsed to d. temlne spectral contrast for any selected wavelr,»g h re- 
solution o; hand Width The J;; 

r',,Vb''u':p:v'\-.p..;;.:i t..o/fpr pu..,.,,., ,.run.ft .a.,.- 


spectral remote sensing missions. A combination of a few selected measurements of solar 
Irradiance and reflectivities, and the subsequent computer calculations of spectral eontrast^ 
could be used to select film/ filter combinations that would yield the highest contrast. 

This method would eliminate costly flights and subsequent data reduction performed in hopes 
of empirically determining optimum bands. 

Calculations From LANDSAT Data 


Modeling of spectral signatures of mixture .- Two approaches for modeling were used in 
this study. The first is a widely reported method that describes a mixture class signature 
in terms of weighted combinations of component class mean vectors and covariance matrices 
(ref. 1, 2, 8, and 13). The model describes the mean spectral response from a ground re- 
solution element as: 




PjMi 


(U) 


'^^re Mp^ ^ response vector for a mixture of N component classes; 

P = Relative amount (proportion) of class i; 

= Mean response vector for the ith component class. 

The above relationship assumes statistical independence of normally distributed data 
points belonging to class i. 

Assuming statistical independence for variables associated with elements from different 
object classes, the relationship for a mixture class covariance matrix can be written as: 


S-f 


(5) 


where Cp = Covariance matrix for a mixture of N component classes; 

P = Relative amount (proportion) of class i; 

= Covariance matrix describing the distribution of the ith component class. 

Equations (U) and (5) represent the model used for automatic classification of mixtures 
of classes with a maximum likelihood processing algorithm. For supervised learning recogni- 
tion processors^ component class mean vectors and covariance matrices must be determined and 
proportions for each possible mixture must be specified. The approach presented here is 
based on modeling the spectral response within a single pixel and should be applicable to 
recognition processing with the RECOG (ref. 10, lU, 15) maximum likelihood algorithm (CLIKE) 
since classification is point by point. 

Most of the work reported in the literature utilizing this modeling technique required 
sophisticated algorithms that calculate various mixture class signatures and select the one 
that gives the closest approximation of the response from a given pixel. In our study, 
however, we used the model to obtain a set of spectral signatures for specified mixtures of 
component terrain classes that are used by the RECOG pattern recognition routines to identify 
all ^6S data points that exhibit that response. The value of this method lies in its straight- 
forward applicability to the existing RECOG processing sequence, thus eliminating the need 
to develop new processing algorithms. 

Equation U was used to simulate the spectral response expected from a mixture of the 
two component classes (A and B) whose extracted mean vectors were shown in Figure 2. Figure 
8 shows the resulting simulated mean vector (d) for a 50^6 mixture of the two classes le. 
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■ P_ • 0,5. A llKoly 0 X|U!«imtion for the difference between curves C end D la the dlffl- 
Ity of eccuretely estimetin^; cover densities from high altitude aerial photographs uaed 
ground control (in this case Mission NASA RB-57 photographs with an approximate 
ale of 1:100|000) and the difficulty of locating non- homogeneous (mixture) training seta 
I gray maps (ref. l6). 

When a set of representative component terrain cJass signatures arc obtained then the 
Ignature of any mixtvure of these component classes ma>’ be determined vising Equations U and 
, A computer program MIX was written to take twc» component class signatures and two propor- 
Lon vectors and calculate the spectral response of any combination. Tlie method involves 
callng the mean vector and covariance matrix of a component class by a proportion factor 
nd adding the result to the scaled mean vector and covariance matrix of another component 
lass. The proport ion factors p; must be within the set 0 < P < 1 and P,i - 1, The re- 
ulting spectraQ signature may then be useil to classify all~plxels in a set of LANDSAT hCS 
ata that exhibit a dnilar spect ral res{x)nse. 

The second modeling technique use^t linear regression to pmlict the spectral response 
»f a mixture- containing pixel leased on the mean response vectors for known mixtures of 
.errain classes. The regression model for a two- component mixture takes the form: 

Yj • C ♦ B(X) (6) 

where Y ■ estimated mixture response in wavelength Mnd 1; 

C ■ constant (Y intercept); 

B • coefficient (slope of regression line); 

X • proportion of one component of the mixture in the pixel scene. 

Mt'ad (ref. 5) Indlcatetl that this method may be useful in estimating the spectral re- 
sponse for varying densities of ponderosa pine, but also notnl that the mtHin spectral re- 
sponse may he affectevl hy the arrangement or distribution of ^wnderosa pine within a train- 
ing set as describe^! by the standard deviation. 


The methods presented helow, with the exception of signatun* extraction, were developt'd 
for mixtures of two component classes. The following Is a synopsis of the techniques used 
to obtain component- and mixtvire-class spectral signatures for automatic analysis of MSS 
ilata : 


METHOD 

TECHNIQUE 

CONltn'KR PRtXiKAM 

Signature 

detraction 

Statistical sam^vllng 
of MSS data points 

KKCOO (nvase Appendix A, P) 

(see ref. lU 15) 

Oemponent Signature 
Estimation 

S<>lving simviltAnt'ovis 
equal Ions 

SIOCALC 


Linear regression 

STATtt'H 

Mixture Signature 
Simulation 

Addition of welghteii 
mean vectors and 
covariance matrices 

MIX 


Signature Extraction .- Wlven areas on the ground have been satisfactorily identified as 
containing a known terrain class or mixture of terrain classes and the graymap coorldlnates 
have been determineii^ then spectral signatures may be obtained by statistically sampling 
these points. This process is known as signature extraction ami represents the conventional 
mode for obtaining signatures that are assumeii t»> characterixe a class of objects. The de- 
signateii data points km^wn to wntain a terrain edass are sampled and the overall mean Is 
determlneti for each v'hannel of Kin vlata. Hvcse mi«ans, four in th»» esse of lANBBAT-1, are 
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collectively Known ns the mean vector. Because of the inherent variability found in natural 
o^cts an additional set of statistics is obtained that characterises the population dis- 
tribution of the sampled data points. This set of statistics is the class covariance matrix 
and describes the within- channel and between- channels variation of the data. The technique 
used in this work is contained in the RECOG routines, specifically Phase 2 (Appendix A). 

The foraulas used in these routines are described in Appendix B. 

Estimati ng component-class spectral sicnatures .- One method for estimating component 
class signatures from known mixture class signatures involves solving two equations with 
two unknowns. The spectral signatures for unknown mixture classes are assumed to be related 
and can be described in equation form in the following manner: 


P^(MSRji 


Pl,(MSR^)l 


^tp, 


(7) 


where P , P * Proportions of classes A and P in the training set, respectively, 


in 


(MSR ) ,(MSIt ). = 

wavelei^th band i. 


Mean spectral responses of the comi>onent classes A and B, respectively 


Mp - Mean spectml response reconied at the scanner for the mixture of classes A and 
B in wavelength band i. 

The terms P . P . and M are Known for each training set^ltviviiig, ttie mean spectr^ re- 
sponses of the eSmpoKent claSses A and B to be detemined. If there f 

spectral signatures describing training sets with similar components but of different proper 
tions, then it is possible to calculate the two component classes by solving the two equa- 
tions simiQtaneously. A general solution for the two-c-omponent class case takes the form. 


^^\i - V/a 


M 


li 




(P) 


a,^ ^1 


- K ^ 


Equation 8 prcluoes the mean spectral response for class P. To Jerlve the mean spectral 
respOMe for class A, the residt can he suhstitut,Kl Into «|uatlon 7 which can he rewritten 

as 


MSK. 


MP. 






If the elements of the covariance matrix for a terrain class spectral signature behave 
in a similar fashion, then this method coidd be used to obtain a calculated t-ovarlance i^ix 
for a c-omponent terrain class. The covariance matrix for a four-channel case such as UNTV 
SAT-1 can be treated as a 10 element army to simplify the calculations. TTiis can be done 
since the off-diagonal elements of the covariance matrix are mirror images of each other. 

A computer program called SIGCALC was written to take two extmeteii mixture 
class spectral signatures and calculate the mimn vector and covariance imtrU 
ponent classes. In addition, the proportions of the component classes in the 
determined from ground truth must be specified. Situations where P. Pp ” ^b^^a, 
avoided since the denominator of equation 8 cannot equal tero. ^ 

The above method for deterrdning component class spectral sigtiaturos should ^ 

liable results if the proportion estimates for the component classes are accurate. 
it is often difficult to measure the proportions of component classes with the 
needed to estimate representative componejit class signatures. Also, 
of spectral responses for component classes found in .uiture it was 
method that estimates the component -cl a.m spectral signatures from mo 


class signatures. This method Involves using a stepwise linear regression analysis. 

The stepwise linear regression analysis used is an applications program available at 
Colorado State University called 3TAT38R (ref. lU, 15). This program was used to develop 
a regression model that estimates the spectral response of a terrain class, either a com- 
ponent or mixture, given a set of mixture-class spectral signatures and the proportions of 
the component classes. 

Simulations of LANDSAT data .- In order to examine the applicability of modeling mixture - 
class sp»*'’^** *l R<gna.t:iirns for automatic classification of LAMDSAT-1 MSS data, a series of 
experiments was performed with simulated data. They include: identifying the expected 

improv«nent of avitomatic classification using modeled mixture terrain-class signatures versus 
conventional component class analysis, comparing modeled and extracted mixture-class mean 
vectors and covariance matrices, and analyzing the signature calculation techni>^ue. 

Creation of simulated data fields: The automatic classification procedure conducted 

for these experiments approximates that of conventional RECOG LANDSAT-1 MSS analysis (ref. 

10, lU, 15) with the exception of using modeled mixture class signatures and simulated 
LANDSAT-1 data. The steps are: l) EIxamlne control data on site to establish locations of 

representative training sets, determining the line and coliunn numbers of the training sets 
on graymaps and extracting the mean vectors and covariance matrices from the data. 2) Simu- 
late mixture signatures with the component class signatures obtained in Step 1 using Program 
MIX. MIX requires the mean vectors and covariance matrices of the component classes as in- 
put, as well as the mixture proportions desired for each mixture class. The simulated mean 
vectors and covariance matrices are punched on computer cards by MIX for future use, 3) Gener' 
ate a random-normal data field for each class in REX^OG fonno.t with program DTAFILE. This 
program uses a random-number generator that selects points as belonging to a given class 
based on the mean spectral response in each LANDSAT-1 MSS band and the appropriate covariance 
matrix, all within a norrAl Gaussian distribution. The overlying assumption here is that 
wildland terrain classes are normally distributed, which conforms to the assumptions imposed 
on the mixture modeling method (ref. 8) and the maximum likelihood pattern recognition algo- 
rithym (GLDCE) of REC(X1 (ref. IC). DTAFILE was designeti to create simulated data fields of 
specified site for each class described by a mean vector and covariance matrix. The gener- 
ated data points are written on a perranent file and stored for later use by any' routine 
t'lat uses data in a RECOG format. 

The simulated data fields for each of the component and mixture classes were generated 
to contain 1000 spectral response values in each wavelength band^with each point being de- 
scribed by four variables (spectral response in each wavelength band). 

The data t'lels as generated by DTAFILE represent LANDSAT-1 MSS data with the added ad- 
vantage of "absolute" ground truth Information. With these fields it is possible to analyze 
the accuracy of automatic classification using the conventional HECOG maximum likelihood 
routine and the mixture terrain-class modeling method. An example data field is shown in 
figure 9 as a graymap of a nine-class fl^'ld. class fields can be seen as horizontal 

bands 100 points across and 10 jHDints down. 

aasslfication of component classes: The initial experiment performed involved using 

a classification of simulated data that contained fields of two component classes and three 
mixture classes modeled from these c*omponents. 'Drfo component classes. Grassland (denoted by 
the symbol, G) and Forest (F), were identified on NASA RB-57 (scale 1:10C,000) color IR 
aerial photographs and located on a griy’map of August 20, 1972 LANDSAT-1 MSS (prame no. 102 - 
17135 ) data over south- central Colorado. Signatures were extmeted and then ust»d to obtain 
mixture class signatures with program MIX. The mixture classes were simulated with a 
fied proportion increment of 0.25, rxx*duclng three mixture classes of 75% Grossland - 25% 
Forest (C), 50% Grassland - 50% Forest (D), and 25% Grassland - 75% Forest (E). The mean 
8 pectral** response curves fer the five classes are shown in figure 10, 


The simulated data fields were then treatcvl as actual LANDSAT-1 MSS data for purposes 
of classification. Component- cl ass signatures were oxtracteii from those fields created with 
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component- class signatures (top 3 rows, fig. 9) and used to classify the entire set of simu- 
lated data. The classification display Is shown in figure 11 with the Grassland field at 
the top and the Forest field at the bottom. The mixture class fields are located between 
the components with C nearest the Grassland field, D in the middle and E nearest the Forest 
field. Note the number of unclassified points (seen as blanks in the display) in each of 
the mixture-class fields. These points were thresholded out at the 11.100 percentile as 
not belonging to either of the component classes. Misclassifi cat ions of the mixture-class 
fields occur only in fields C and E, with points being miscla&sified as the component class 
which has the highest proportion in the mixture. 

The absolute classification accuracies were determined for each clas . Since each data 
point was created as a specified class, a point was considered misclassif ied if it was classi 
fied as any class other than that specified. The results are summarized in the following 
classification confusion matrix (CCM) listing each class, the number of points that were 
classified correctly, and the number of points misclassif ied as other classes. 

CLASS 

FOREST (F) 

2 % GRASSLAND - 
FOREST (E) 

•jOi GRASSLAND - 

50i FOREST (D) 

GRASSUND - 

25% FOREST (C) 

GRASSLAND (G) 

The table is read with the true class name and symbol on the left and the symbols of classes 
that a point was classified as across the top. The diagonal elements of the matrix represent 
accurate classifications. The off-diagonal elements represent either Type I or Type II errors, 
defined as follows. 


F E D C G 

1000 

1U8 

376 

1000 


Decision 

Classify AS A 
Do Not Classify AS A 


True Classification 

X belongs to A X Does Not Belong to A 
Correct Decision Type I Error 

Type II Error Correct Decision 


The overall accuracy of this classification was only 40.0% which in most LANDSAT-1 appli- 
cations should be considered quite poor even though the component classes were classified 
perfectly . 

Classification of two component and three mixture classes: The objectives of the study 

presented in this section were twofold. The first was to verify that the mixture- class 
spectral-signature modeling technique produces signatures that can be used by RECOG to classi- 
fy pixels that contain mixtures. The second, assuming the first to be satisfied, was to de- 
termine the increase in information acquired using modeled mixture-class signatures over 
classification conducted with only component- class signatures. 

The same data-field generated for the analysis of component classes was used for this 
experiment. Treating the data file as actual LANDSAT-1 MSS data, component- class mean vectors 
and covariance matrices were extracted. Signatures for the mixture-class fields were not 



extracted from the data, but Instead vore ir.ovicled from the component- class mean vectors and 
covariance matrices. Hiis procedure follows that which could be used with actual IAKDSAT-1 
data since representative mixture- class training sets are usually difficult to locate in 
natural situations. The two extracted component and three modeled mixture- class signatures 
were used to classify the simixlated data file. 

At this point it was desirable to compare modeled and extracted mean vectors and covari- 
ance matrices. The signatures cf the mixture classes were extracted from the simulated data 
and canpared with those that were modeled from the two component classes, as follows; 
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BAND 

MEAN 

STD. DEV. 

CUSS 

BAND 

MEAN 

STD. DEV. 

CUSS 

BAND 

MEAN 

STD. DEV. 

CUSS 
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CUSS 
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The results indicate that the cxtracto<i and modeled mean and standard deviation vectors vary, 
at most, by less than 0.1 of one standani deviation unit, which can be considered as quite 
adequate for this analysis. Figure 1C shows the classification display of the results. A 
classification confusion matrix (CCM) was also civnpllod for this analysis. It Is a comparison 
of extracted and modeled mean and standani deviation errors. Tlie standard deviations shown 
represent the square root of diagonal co\ririance elements. The results are as fcllcws; 


CLASS F 

FOREST (F) 9^ 

C5% GRASSLAND - q 

FOREST (F.) 

50^ GRASSLAIiD - 
yyfi- FOREST (D) 

75i CRASGUND - 
FORFirr (C) 


E 

3 

<rn 

8 


D 

C 

G 

14 

o(>7 

25 


12 

vr*o 

18 


Thr classification results Indicate that simulated mixture-class spectral signatures 
can be used to classify MSS data using the RECOG program. The classification information 
and accuracy were also greatly Increased by the use of mixture signatures. The overall 
classification accuracy using this technique Increased to over with only a slight de- 
crease in the accuracies of component classes. It should be noted here that automatic classi- 
fication of generated data fields will normally produce high classification accuracies due 
to the lack of "alien objects" which are often found In real life situations (ref. 3)* The 
accuracies produced here should^ instead, bo considered a measure of the effectiveness of 
the classification technique. These results do Indicate that given a pixel containing a 
mixture, It should be possible to accurately Identify It with this method. 

Classification of three component and six mixture classes: The successful results ob- 

tained with the five-class analysis warranted further study, so an additional component class 
was added to the signature set. A mean vector and covariance matrix for Water (W) was ex- 
tracted from August 20, 1972 LANDSAT-1 MSS data. Mixture signatures were modeled from the 
Grassland component class and the new Water class with a proportion Increment of 0.25. The 
resulting three mixture classes and the component class combined with the original signature 
set bring the total number of signatures used to nine. Figure 13 shows the mean spectral 
response curves for the Grassland-Water mixture classes. The separation of these classes 
was adequate in all bands. The mean spectral response curves for all nine classes are Illus- 
trated in figure l4. 

The classification procedure followed that of all previous experiments with a new data 
field being generated to Include all nine classes. The classification display Is shown In 
figure 15. A classification confusion matrix (CCM) for the 9-clnss analysis Is as follows: 
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50% WATER (R) 
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25% GRASSLAND - 
75% WATER (S) 
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WATER (W) 






962 


These results show a good classification accuracy for all the classes, with an average classi- 
fication accuracy of 98^t. Note that the greatest confusion exists between mixture classes 
that have the same proportion factor for a given component class, l.e., 75^ Grassland-25^ 
Forest and 75lt 0rassland-25lt Water. This Is probably the result of masking of the lesser 
components by the Grassland response. 


Thr* following comparison of th# RKxlolod and extracted moan and standant dovlations was 
mad* as an additional check of the mixture modeling method: 
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As In the iNise of the five-class classification! the simulated mixt.ure signatures closely 
approximate those extracted from the generated datafiel.ls. Ttu'se results imply tl»at if a 
representative set of com}v»nent signatures can he found, then mixtures of these classes may 
he successfully mi'deled anil used to classify pixels containing these mixtures. 

Test of signature estimation technluuo! In wildland ar«N\s where large arivin of homogane- 
ous terrain classes are dlffleuli to locate, it becomes ni'cessary to id'tain representative 
spectral signatures from, a relatively small numher of pixels. I’u some areas it is even dlffi- 


cult to find any pixels that contain only one terrain class (ref. l6). Consequently, methods 
to obtain component terrain-class signatures from pixels containing mixtures of terrain 
classes In Known proportions have been Investigated (for example, see ref. 17). The method 
involves solving a set of llncur equations assuming that the data is from a Gaussian distri- 
bution and each signature has a common covariance matrix. Our study used a similar, but 
much simplified technique to attempt to acquire component-class spectral signatures from 
training-set data of terrain classes of Known mixtures. A test of the method to obtain 
component- class spectral signatures using equations 8 and 9 was conducted using slmulateti 
data fields created in the same manner as the preceding experiments. The specific objectives 
of this test were; l) to verify that this method could be used to determine component class 
signatures suitable for automatic mapping of LANDSAT-1 type MSS data when only mixture ter- 
rain classes were available for signature extraction, and 2) find out the type of classi- 
fication accuracies that can be expected when implementing these estimated signatures as 
input into the model and RKCOG. 

The initial test used the random-normal data fields for the nine-class analysis described 
above. Spectral signatures were extracted from the data for mixture classes 75^ Grmssland- 
^ Forest (C), ^ Grassland- 501^ Forest (D), 7 ^% Grassland-25% Water (q), and 50% Grassland- 
50% Water (R). The pairs of extracted signatures and their corresponding proportion factors 
were input into SIGCALC and the component signatures for Grassland, Forest and Water were 
calculated. 

Comparison between component class mean and standard deviation vectors extracted from 
the data set and those calculated with SIGCALC is as follows: 
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The calculated component class signatures closely approximate the extracted^ with the average 
deviation being 0.08 for mi*an vectors and 0.15 for standard deviations. 

In order to compare the effectiveness of estlinated signatures with extracted signatures, 
the nine-class data field was classified by the computer. 'Hie class signatures used as in- 
put included the three estimated component signatures and six mixture-class signatures model- 
ed from the estimated signatures. The classification accuracy obtained was approximately 2% 
less than that obtained when using extracted component signatures. Hie classification map 


^ is shown* in figure l6. 


follows : 


The classification confusion matrix (CCM) of this data set is as 
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An overall classification accuracy of 96.1^^ was obtained, indicating that good classification 
results can be expected using this method. ciassiiication 


The simulated data fields that contained the mixture-class signatures used to calculate 
the co^nent classes represent uniform data. That is, each data point (pixel) represents 

objects. Because this situation is sometimes difficult to find in natural 
oi signatures extracted from non-uniform data 

of Non-uniform training sets can be defined as areas that contain a mixture 

If the^orL^fo^^t single pixels containing various proportions of the components. 

If the non-uniform training set is considered in its entirety the mean spectral resnonse 

training set with the same proportions of component 
?0 uniform and non-uniform data is illustrated in figure 17 as two 

Sid^f ® ® unifonn mixture of two component classes (white and 

1 containing a non-uniform mixture (b), but both with identical overall 
blaS ^ component^asses. In the uniform training set each pixel (outlined by heavy 

theintlrrL^iniM^Lt component classes, producl!;g a 50 % mixture with iS! 

Rr»nn«r Spectral response from a pixel is integrated at the 

^)^h^"“cell distribution of the components does not affect the between- 
coitai^liri^2^ training set. The non-unifom training set is composed of some pixels 

containing 100^ of the other class and some that contain 
various mi^ures of the two components. The overall mixture proportion averaged over the 
entire training set is ^ of each component class. Referring to the equation 28 (Appendix) 
easily deduced that the mean vectors obtained from each training sot should equal ^ 
each other assuming constant illumination and no within-class variabUity. By equation 30 
^Appendix) however, It can be seen that covariance matrices will most likely vary due to the 
betveen-pixel variability existing in the non-uniform training set. 


To verier these statements^ mixture-class signatures from simulated non-uniform data were 




signatures. A comparison of extracted and calculated mean and standard deviation 'vectors of 
the component classes for uniform and non-uniform(*) data are as follows; 
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Note the close agreement of means and standard deviations from the uniform data and the close 
agreement of means but non- agreement of standard deviations calculated from the non-uniform 
data. 


The consequence of using such large covariance matrices obtained from the non-uniform 
data is illustrated in Figure l8. The nine-class simulated data fields were classified with 
the calculated component signatures listed above. The first classification map (fig. l8a) 
was obtained using signatures calculated from uniform training sets and the second (fig. l8b) 
was obtained using signatures calculated from non-uniform training sets. The classification 
accuracy of component classes obtained for the uniform case was 99*1^ whereas the accuracy 
obtained for the same classes in the non-uniform case was only 33 The latter figure is 
misleading in that the nine-class data set was classified as 97*9^ Grassland, 1.7^ and .4^ 
water. These figures indicate that the class distribution for Grassland was improperly de- 
fined by the covariance matrix. 


The above study, coupled with the frequent lack of uniform mixture-class training sets 
in natural situations, lends itself to implementing an artificial covariance matrix for de- 
scribing distributions of component classes. The use of a covariance matrix common to each 
class was discussed (ref. 8)for another, more sophisticated, signature calculation technique. 
Mean vectors for component classes can be successfully calculated from non-uniform mixture- 
class training seta. If a common covariance matrix that describes the proper distribution 
can be found, then with the calculated mean vectors, the signature should be representative 
of its respective class. To verify this, the component- class mean vectors calculated from 
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non-unlfom mixture class training sets and various common covariance matrices were used to 
classify the simulated nine-class data fields. 

The first classification utilized a covariance matrix with diagonal elements ( ff.. ) of 
(1.3i3.l) and off-dlagonal elements of 0.5. 'Hie determinations for the covariance elements 
were made by a priori Inspection of covariance matrices for component and mixture classes 
extracted from simulated uniform data. The standard deviations In each band represent an 
estimate of the average standard deviations for all classes. The off-diagonal elements were 
established to be 0.5 as an estimated average of the correlation values for all of the 9- 
class correlation matrices. 


The classification procedure followed that of all earlier runs with the mixture class 
mean vectors being modeled from the component class mean vectors calculated from non-uniform 
data. The classification map is shown in figure 19; its corresponding COM Is as follows: 
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17 

•?59 

15 


9 

2 

25% GRASSUND - 
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21 
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75% GRASSLAND - 
25% WATER (Q) 



20 

23 


956 



50^ GRASSLAND - 
50^ WATER (R) 


25^ GRASSLAND - 
75i WATER (S) 


1000 


WATER (W) 


1000 


The acceptable level of classification accuracy cbtained (‘J7.3^') indicated that utilizing a 
common covariance matrix iray help solve the problem of establishing representative spectral 
signatures for canponent classes when there exist Insufficient points for conventional sig- 
nature extraction. 


An additional classification was conducted with the same class mean vectors as the above 
analysis, but with a common covariance matrix with the same diagonal elements and off-diagonal 
elements set to 0.0. The purpose of this run was to establish the expected classlfl^tlon 
accuracy for 0.0 correlations between channels. The classification results were again pre- 
mlslng^wlth the average accuracy obtained being 96.5%. The classification map is shown in 
fig. 20; its corresponding CCM is as follows: 
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23% GRASSLAND - 
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75% FOREST (C) 
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25 

66 
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983 


73% GRASSLAND - 
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22 
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30^ GRASSIAND - 

WATER (R) ICXX) 


23^ GRASSLAND - 

1% WATER (S) 1000 


WATER (W) 


1000 


One additional common covariance matrix was calculated, and tested as suggested by 
Nalepka and others (ref, 8, p, 17). This matrix consisted of an average of the covariance 
matrices of their signature set. One immediately recognizable disadvantage of this method 
-s that a set of representative covariance matrices must be obtained from the MSS data which 
have small values such as those expected from component class covariance matrices extracted 
from uniform data. The results of using this technique for calculating the diagonal elements, 
but setting the off-diagonal elements to 0.0 are shown in figure 21. The corresponding CCM 
is as follows: 
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Note that the 9- class data classification results are slightly less {9^. Off,) than those ob- 
tained using an estimated average of the diagonal elements and setting the off-diagonal ele- 
ments to 0.0. 

The overall results obtained from using common covariances appear promising for classi- 
fying LANDSAT-1 type MSS data when signatures must be obtained from non-uniform mixture class 
training sets. This technique applied to automatic classification of actual wildland data, 
is described in the next chapter. 

Tests using LANDSAT-1 MSS data .- E>cperiments similar to those performed with simulated 
MSS data were conducted with actual LANDSAT-1 MSS data. An initial experiment was performed 
using extracted assumed "pure" component-class spectral signatures and mixture-class spectral 
signatures simulated from these component classes to classify an area designated as the 11- 
Mile study area. An additional experiment using estimated component signatures and simulated 
mixture signatures was performed using several plots located in the Manitou study area where 
ground control data were more detailed. Both of these experiments utilize methods described 
above. The objective of these experiments was to relate the concepts Just presented to a 
real-world application. Hopefully, this will provide insight as to how well these techniques 
work to improve classification accuracies and/or what other factors may need to be considered 
to obtain optimum results. 

LANDSAT-1 MSS data acquisition and description: The set of LANDSAT-1 MSS data analyzed 

is from an August 15, 1973 (frame number 1388-1713^) overpass of the test site. The data 
set was provided by Dr. Richard S. Driscoll, Remote Sensing Project Leader at the U.S. Forest 
Services Rocky Mountain Forest and Range Experiment Station located in Fort Collins, Colorado. 
The data were geometrically corrected (deskewed) by the Laboratory for Applications of Remote 
Sensing (LARS) located at Purdue University, West Lafayette, Indiana (ref. 5, l8). The ad- 
vantage of the deskewed data is its better conformity to aerial photographs, especially when 
displayed as a microfilm graymap or recognition map. A large section of this LANDSAT-1 frame 
was also studied by Heller and others (ref. U) and Fleming (ref. 19). 

Two areas from this frame were selected for study. They are referred to below as the 
Elevenmile Canyon Reservoir Study Area and the Manitou Study Area. 

Experiments were conducted in these studj' areas using extracted component class and 
simulated mixture-class spectral signatures to classify a portion of LANDSAT-1 MSS data. 

The steps involved in the classifications are; (l) selecting component classes, (2) extract- 
ing representative spectral signatures, (3) simulating specified mixture-class spectral sig- 
natures, (4) classif^-ing the data with a maximum likelihood processor, and (5) analyzing the 
results . 

Elevenmile Canyon Reservoir Study Area; This area was .elected because it is represent- 
ative of the diversity of topography, geology, and vegetative cover throughout much of the 
total test site; and because of the availability of LANDSAT-1 MSS tape data. The location of 
the site is shown in figure 1. Elevations range from 2500 to 3020 m. The site is at the 
southwest margin of South Park: it includes the eastern part of Elevenmile Canyon Reservoir. 

Mountainous terrain lies to the north and cast; a wide rolling grassland area lies primarily 
to the south of the reservoir. Wet meadows occur primarily along stream courses. 

Two principle types of rocks occur; volcanic rock and granodiorites. The most commonly 
found volcanics are upper and lower andesite members of the Thirty-nine Mile volcanic field. 
The upper mend>er is generally associated with high flat-topped mountains and the lower with 
the low rolling hills characteristic of South P&rk. 

The granodiorite found in the site is dark gray, medium to coarsely crystalline and con- 
tains gneiss. It is correlated with the precanbrian Boulder Creek Granodiorite of the central 
front range of the Rocky Mountains. It is found mainly in the northern portion of the test 
area where it forms high steep mountains. Several outcrops can be found to the east of Eleven* 
mile Canyon Reservoir and on both sides of Elevenmile Canyon, but their aerial extent is limit- 
ed. Elsewhere, most of the bedrock is covered by some form of vegetation. 


<n»e vegetation In the area consists of a variety of forests and grasslands. The forests^ 
which cover approximately one-third of the test site^ Include; (l) ponderosa pine, (2) Douglas- 
fir, (3) spruce fir, subalpine fir, and (4) Aspen. According to Heller and others (ref. U), 
these forest species occur as "pare" types, but in many places there are various mixtures of 
species. There also exists a wide range of tree canopy densities, from very open to so dense 
that crowns nearly touch. In open stands of trees enough sunlight exists to allow the de- 
velopment of extensive herbaceous and/or shrubby understory, but under conditions of dense 
tree cover, little understory cover has developed. 

The principle understory vegetation associated with Ponderosa pine forest includes 
Arizona fescue and mountain muhly. At higher elevations these species give way to Idaho 
fescue and Thurber fescue and oatgrass. 

Approximately one- third of the study site consists of low rolling hills covered with 
vegetation characteristic of the short-grass prairie found throughout South Pirk. The prin- 
ciple species are blue grama and slim stem muhly. Mountain bunchgrass communities can also 
be found at the interface between grassland and forest. 

Also occuring in the area are wet meadow and stream bank communities. These communities 
occur along the shore of Elevenmile Canyon Reservoir and along streams. Various species of 
sedges, rushes, and bulrush occur either as monospecific or mixed stands in the moist areas. 

In drier areas, blue grass and tufted halrgrass are found in mixed communities. The promi- 
nent shrubby communities consists of willow and shrubby cinQuefoil. Although they occur 
sparsely throughout the study area, these communities are generally found associated with 
the wet meadows. 

Five component classes were selected for use in this study area. They are Forest, Grass- 
land, Water, Mountain Grassland and Wet Meadow. These classes were selected because of their 
importance in the area and also because mixtures of these various types occur. The Forest 
class is composed of several types of trees with a minimal amount of understory vegetation 
showing through the canopy. The Grassland class is characteristic of the short grass prairie 
community found throughout South Purk and usually has much bare soil associated with it. The 
Water class is representative of the water of the reservoir. The Mountain Grassland class 
represents grassland communities found at higher elevations and at the Interface between Forest 
and Grassland. The Wet Meadow class was selected to represent irrigated meadow areas and 
stream bank communities. 

Homogeneous areas representative of each class were delineated on aerial photographs. 

These areas (training sets) were then transferred to LANDSAT graymaps and the line and point 
coordinates were detenr.ined. Figure 22 shows the microfilm graymap of MSS BAND 5 for the 
study site. Some modification of the photo locations were made on the graymaps as dictated 
by the graylevels associated with each training set. In this vay^ anomalous data points were 
avoided before signatures were extracted. A number of training sets were located for each 
class to provide statistically representative signatures. T»'e number of points used to cal- 
culate each class signature were; Forest (75), Grassland (26), Mountain Grassland (20), 

Wet Meadow (20), and Water (96). The number of points for Grassland, Mountain Grassland and 
Wet Meadow was small due to difficulties encountered In locating training sets on computer 
graymaps, but are assumed representative for this study. 

Mean vectors and covariance matrices were extracted for each training set for each class. 
The mean vectors for a class were compared with one another and the eigenvalues and eigenvec- 
tors of each class training set were determined and analyzed. Those training sets exhibiting 
dissimilar mean vectors and eigenvector plots were discarded so that final signatures obtain- 
ed were ns representative as possible. 

The final component- class signatures were then used to classify the area. The resulting 
classification map is shown In Figure 23 ns n microfilm display. As a check of how well the 
signatures Identify their respective classes, the number of correct classifications In each 
training set was determlnrd. All of the clnsses had training-set classification accuracies 

of 100% except Water which had 99%» A total of 60.5% of the test area %rae classified as one 

of the indlMtiiic that for a oonventlonal analysis such as this, more elasaes 


To antisfy this need, the following 


2 . 


3. 


nee»l to bo added in onier to classify the entire aro«i. 
mixture classes were selected; 

1. Grassland-Forest select»*«l to identify pixels covering areas with trees growit\g 
over Grassland areas and interfaces between the two cownunlties. 

Grassland-Wet Mt^idow-- selected to identify pixels falling on the interface of the 
grass covered areas and Irrigateti meadows and stream bank c'ommunitles. 

Forest- Mountain Grassland--aelected to identify pixels covering the Interfaces be- 
tween the two communities and locate areas of Mountain Grassland that appear as 

openings in the Forest which are smaller ♦■•han the resolution of the LANHSAT-l MSS 
scanner, 

A proportion increment of 0.25 was selected for each combitjation of classes, and mixtureJ 
class spectral signatures were simulateii with program MIX. The final signature sot useii to ^ 
classify the teat area included the five extracted comjx^nent classes and nine interme^ilate 
mixture classes. The micrcfilm classification map at the 11.100 Chi-Square threshold level 
is shown in Figure 24, An incr«\sse in the number of classified points resulted in 86 
of the ^ea being cUssifi.M, mthrr than Mi.st of the unclassified points are kssooi- 

atod with the interface between Grassland and Water since no spectral signatures character- , 
istic of these mixtures were used in the classlfiivition. Classification of Grassland-Wet 
Mea^ and Ptorcst-Mountatn Grassland mixtures ap(>car rtvisonablc. Grassland-i'\?re8t mixtures 

Grnasland-?5'S> r\>rest also appi'ar to be classified well. 

The 75^ 0’*r.';slsnd-c5^ Ft^rcst mlxturt' class i'ontained a lai-ge tum'unt of obvious misclasslfl- 
catlons in artvis where no Forest was present. 

To ^Investigate the source of these mlaclassifications, spectral signatures were extract- 
ed from the data for ar.'ns classified as the 7^^ Grassland-251t. b^orest class. A comparison 
or the mt*an and standanl deviation vectors for t»>e cxtractiHi signature ajul the simulated mix- 
ture slgkiature showcti then to be almost identical, as shown below; 


meted 

MSG RAND 1234 
Mean 20.4.' 28.? > 37.34 19.21 

Std. D«*v. .‘)li l.«xj 1.51 1,00 


Simulated 


1 

28.0? 

l.tV* 


2 

28.87 

2.02 


3 

37.19 


4 

IQ. 00 

1.15 


indimtes that the ratorlal tl.at is being classlfl.M as a Grassland-PV'rcst Mixture 
has a vegetation or soil ty]M' that prv>ducea a resptMisr similar to tliat calculated for the 
e ^ronsequenees of such a mlsdassl fKvxtlon are an overesliration of the an.ount 

of Forest. In the scene and an underestimation of Grassland, because of the difficulties cn- 
ccuntertHi in establishing mixture den.sltles frvvr. high .altitude aerial photograplis no quanti- 
tative accuracy analysis was obtaimvl for tt»ls experiment. ^ 

IW'cause of the gro.-.s nlsclasslflc.atlons of the 7% Crassland-2'-4 bVrest class, the nrt>- 
portton Increment used to .-^’miaate mixture clar.s signatures was char;gesi to 0.33 and iu*w sle- 
jMitures were slmulativl. The residtlng signature set included the five co«nponent classes and 
ten two-comix'nent mixture classes, combinations of Grassland and Water and 'Water and Wet 

Meadow were made as an attempt to idenlifi- those pixels left, unclassifl.xi aianind the reservoir 
in the previous experiment. 

A classification was mde^proiiucitvg the mlcrt'fllm display shown In figure 25. The re- 
sults Indicated that the anas previously mlsclassifled as 7H Grassland-?*!^ ft^rest were still 
i?*^***^***^^^'^ idenUfled. A number of pt^lnts along the shoreline of the reservoir were 
Identiflisi as either o7t Grassland- 3?1f.water or 3JJt Graasland-67% WKter. There appearnl to 
be s«ae confusion of 67^ Wet Meadow-331t Whter and 331t Wet Meadow-671. Fhrest as tbreomponent 
Class with the greatest amount oecnrlng In the mixture. The ether mixture classes appear to 
be classified correctly- since they appearnl where mixtures would logically exist. 

Manltou studj* area; Tt»e Manltou study area Is 
C,v.,n R..,rvolr .r.vi 1„ ,ro,.„,i tho U.S. Sorvlcc mnHo. 
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a fev selected plots identified by Mead (ref. were studied in this work. 

The plots lie within an area located between 38 * 52* 30” and 39* 22* 32 " north latitude 
and 105 * 00" and 105* 37* 30" west longitude, along and just north of the area shown in fig. 

1. The area typifies the mountainous part of the east slope of the Rocky Mountains, with a 
large variation of slopes and aspects. The topography is dominated by mountains and valleys 
with a north-south orientation due to extensive faulting that occured during the Pliocene 
epoch and earlier. 

For the most part, geologic materials in the area are covered with vegetation. The 
principle materials are pinkish to reddish soils developed from Pikes Peak granite. The 
underlying geologic unit is the Fountain Formation, but exposures are limited. 

Variations in topography in the area have resulted in a complex mixing of plant commun- 
ities generally aligned with the topography. For the most part, the plant communities are 
similar to those described for the forested hillslopes, mountains, and wet meadow areas of 
the Elevenmile Canyon Reservoir study area. 

In this study area, experiments were performed to classify LANDSAT-1 MSS data with com- 
ponent class signatures estimated by solving simultaneous equations and linear regressions. 

Both of these techniques require data where the component classes and their respective pro- 
portions of the scene are known. 

Several plots located in the Manitou study area and identified by Mead (ref. 5) were 
analyzed to obtain a number of plots suitable for signature estimation. Mead was specifically 
interested in the percent cover of ponderosa pine of the plots so it was a logical choice for 
one of the component classes. The second component class was identified as the other materials 
in the scene and denoted as Background. The Background class is actually a mixture in itself 
of soil, bare rock, and understory vegetation. Use of this Background class necessitated 
using only those plots wiiich contained similar species occuring in the same relative amounts. 

It was also necessary to identify plots with similau* slopes and aspects to avoid the signature 
variability associated with those factors. In summary’, the criteria used for selecting the 
plots used for signature estimation were: 

1. The plots must have known percent covers of ponderosa pine. (The Background percent 
cover was determineil by subtracting the ponderosa percent cover from lOO). 

2. The plots must contain similar l>ackground materials such as vegetation, rocks and 
bare soil. 

3 . The plots must have similar orientation as determined by slope and aspect. 

The first criterion was met by analysis of aerial photography by five photo interpreters from 
the U.S. Forest Service. Hie second was accomplished by field studies by Mead using a line- 
transect method to establish the materials present in each plot as well as their frequency 
of occurance. The third was met by locating the plots on topographic maps and determining 
slopes and aspects by a computer program called TOPCXlO^ad , ref. 5 * Shiurpiuick and Akin, ref. 
20 ). 


Analysing all of the plots with these criteria produced a small subset of only six plots 
representing 120 total LANDSAT-l MSS data points. The plots are as follows: 
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PLOT 

NUMBER 

MEAD’S PLOT 
NUMBER 
(ref. 5) 

PONDEROSA PINE 
PERCENT COVER 

1 

1 

46 

2 

2 

74 

3 

3 

76 

U 

6 

64 

5 

lU 

22 

6 

17 

16 



SPECIES LIST 


SLOPE 

ASPECT 

(see 

FREQUENTLY 
OCCURING 
BACKGROUND 
VEGETATION 
species list, celow) 

7.9 

96 

J, B, M, P, Q, D 

17.0 

53 

B, J, M, P, D 

6.9 

lo8 

J, B, M, Q, P, D 

5.0 

59 

B, M, J, Q, F, D 

8.1 

77 

B, M, F, P, B 

9.9 

134 

M, AF, P, Q, BB 


ABBREVIATIOU 


NAME 


F fringed sagebrush ( Artemisia frigida Willd) 

B bear berry ( Arctostaptiylos uva-ursi (L. ) Spreng.) 

AT Arizona fescue ( Fostuca arizonica Vasoy) 

J common juniper ( Juniperus communis (L.)) 

P prairie junegrass ( Koelaria crista ta (L.) Pers.) 

M mountain muhly ( Muhlenbergia montana (Nutt.) Hitch.) 

Q quaking aspen ( Populus trenuloidcs Michx.) 

D Douglas-Fir ( Psoudotsuga menziesli var. glauca 

(Beissn.) Franco) 

BB bottlebrush squirrel tall ( Sltanion hystrix (Nutt.) 

J .G.Sm. ) 


Spectral signatures for the six mixture plots were extracted from the LANDSAT-1 MSS data. 

The mean spectral response curves for the six plots are shown in Figure 26. Several pairs 
of signatures along with their respective mixture proportions of ponderosa pine and background 
were used as input for program SIGCALC. The results showed a vide variation in the estimated 
mean vectors for both of the component classes. In addition, the estimated covariance matrices 
contained very unrcGdlstic values including negative variances. Consequently, this method 
%ras found to be unreliable for estimating component class spectral signatures for this case. 
Since results obtained with simulated data described above were acceptable, the error was 
probably due to inaccurate estimates of percent cover and variations of materials associated 
with the Background class. Since no representative spectral signatures could be obtained, 
no modeling of mixture class signatures or classification of MSS data was attempted. 

Since the signature estimation technique used in the above study uses only data for two 
mixture training sets, slight errors in percent cover estimates will result in large errors 
for the calculated signature values. Also, the signatures, if successfully estimated, may 
only be representative of the material found in the specific training areas. To compensate 
for this all of the mixture training sets were used to detennine component class signatures 
by means of linear regression techniques. 

The method used for estimating comptmen'. lass signatures with linear regression was 
discussed above. The necessary input data wore the mean spectral response recorded from each 
of the plots in each wavelength band and the proportions of the two component classes, ponderosa 
pine and Background, tiodc-ls were developed for each wavelength band as follows: 


MSS BAND 

Regression Model 

.5- .6 

MSR=2U.07-5.98(Pp) 

.6-. 7 

MSR=20.9e-7.60(Pp) 

.7-. 8 xm 

MSR*32.05-7.1*U(Pp) 

.8-1.1 -wm 



where MSR= Estimated 
Pp»Proportion 


mean spectral response 
of ponderosa pine in the scene 


No regression model could be developed for the .8-1.1 m wavelength band so a review of 
estimates of percent cover for the plots was made. It was found that for sane plots, the 
cover estimates made by the five photo interpreters varied by as much as 30^. The standard 
deviations for the mean cover estimates ranged firom 5.** to 13.0 for the six plots. Additional 
cover estimates were then made using a dot grid superimposed on enlarged photos of the six 
plots. For these estimates no species determinations were made, but all trees were counted 
as belonging to a Forest class and all non-treo materials were counted as RicKground. New 
regression models were developed as follows: 


yjSS Band 


Regression Model 


.5- .6 
.6-. 7 -wm 
.7-. 8 <*m 
.8-1.1 aim 


MSR=05.70-10.22(Fp) 

K.SR=23.J*6-ll4.03(Fp) 

MSR=3i?.2^-13.23(Fp) 

MSR»19. 55-8.73 (Fp) 


where MSR* Estimated mean spectral response 
Fp* Proportion of Forest in the scene 

High R^ values Indicate that tiie grid estimates for the plots are reasonable values for 
the amounts of Forest and Background. Since only six plots were usetl in developing the re- 
gression models they nay only be applicable for analysis of these plots. Component- class sig- 
nature estimates can be detcrmin«*d from the models by setting the Forest proportion (F^) to 
1.0 and calculating the MSR's in each band to determine the Forest signature and setti^ F_ 
to determine the Background signature. The mean spectral response curves for the estimated 
Forest and Itaekground Signatures are shown along with the extracted mixture curves for each 
of the plots in Fig, 26, 

The estimated mean vectors for the Forest and Background classes were input into program 
MIX and mixture class signatures wore simulated. The proportion increment used was set at 
0.25, producing three intermediate mixture signatures. No covariance matrices for the com- 
ponent-class signatures were estimate*! by the regression models due to the non- uniformity of 
the mixture plots, so a common covariance matrix with diagonal elt*m.ents of 1.0, 3»0» 3*^» 1*0 
and off-dlagonal elements of 0.0 was usoti for each of the classes. 

The six plots selected for the stud^- were classified with a maximum likelihood pattern 
recognition processor. Initially, only the estimated signatures for Forest and Background 
were used to classify the plots. The total number of points classified as these classes was 
only five of 120, or approximately U^, Illustrating the need to account for pixels that con- 
tain mixtures, A second classification was coTiducted with the two component class signatures 
plus the three 8lmulate«i mixture class signatures. The number of classified points Increased 


to lir, or approxlnviloly 9Hi '^Ith Implom nt.ntiou of tho ir.lxiurt -olnss sipnnturos. 


Tht' tolfil xrioant of »'ach I'OK'ponrnt clnas in the plots was dotenniiirii by siunmlng the 
number of points clnssifio4 as a class a!ul multiplying by the propt'rtion of »\ach component 
v'lass. A comparison of the nercent cover estimates for Forest for »'ach plot ns determined 
by the dct-grid method and the classification map are as follows; 


PLOT 


1 


3 

U 

5 

6 


C.HID COV>:S KLTLV-ATKCit) 


55.0 

65.0 
60.0 
UC.O 
00.0 
X ).0 


CLAi,?SIFICATlON 
COVER ESTIMATE 

ii) 

63.75 

71.0*^ 

71.05 

l>f^.75 

0?.75 

37.50 


If the dot-prld cover estimates can be assumt'd accurate then the apx’>unt of hVrest occuring I 
the six plots was over-estimated by ll.'*t. 


EITECT OF THK ATMCiJlIiKKE ON ^^FlXTlvAL aFFLECTANOE 


Kemt'te Sensing I'f TJirget Slgtiatures 


Ttie signature of a target, 
F.RKP), is a radiance signature. 
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N , given bv 
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where H is the total (l^''*, global) incoming solar irratlance, T is the beevr. transmittance 
(either the reflected target beam or the solar bran) of the atmos]'here, is path (between 
the satellite sensor and the target) radiance intrcKiuced b.\ atpjospheric scattering of sun- 
light, and P is the target rtflect ivlty. Tliis target rt'flectlvity, P , is the true signa- 
ture of the target. As can be seen in evjjatlon (lO), the target's radiance signature is 
made up of three fwaroraeters other than the true target reflectivity signature. These para- 
neters are the total incoming solar Irradlance-H, the atuacspiierlc beam transmittance-T, and 
the atmospheric path radlance-Np. Ttuj will be discussed In the following section, where 
all three will be referred to as atra'snherlc effects. Tive vtvelength rt'glon censldert»d in 
this paper is the solar refleetlve region from to l^CO ?vinometers(tun). This region in- 
cludes the lANIVlAT bands (S'C-t'CV, taX'-700, VV'>-t\X'', and ‘VV-1100 nn) and most of the SJkylab 
F.REP S-190 scanner bands. 


B«'oause the target radiance signatun' is depend*»nt upon almostxjeric effects, one goal 
of this program is to I'orrect the satellite uata for these atmosplu'rie effects. vVmparlsons 
can then be made concerning whnt imprv'V''ments can be realired in autvvnatic terrain mapping 
by computer by using algorlthmus for ^x'rrectln*: for attvspherlc effects. 

It is not within the scope of this p.»per to discuss each of these pammett'-*s in great 
detail in terms of atmospheric physics (radiative tran.*sfer); instead, the folliarflng sub- 
sections discuss thtmt In sufficient detail to relate ther to rervte sensing. 


AtnospJierlc !V»air. Trnnsnlttan *e 


At im spheric boom tranaaittaisee, 
Lambcrt's,or Pouguer's l iv) 


T, can t'C given by the folU'Wing expression (BU'ers, 



where Ip Is the path length, relative to the vertical path length of the beam, and T Is the 
atmospheric optical depth. The atmospheric path length is given by sec 6. or sec 6, where ‘ 
e» is the solar zenith angle, sec 6, is the path length of the solar beam, d is the nadir 
angle of observation, and sec G is the path length of the target reflected beam (see fig. 27). 
At low sun angles or oblique observation angles, 9, or G > 70*, the path length needs to be 
corrected for atmospheric refraction (ref. 21). The atmospheric optical depth is wavelength 
- dependent (?v) and made up of three components; molecular (Rayleigh) scattering,'?^ ; part- 
icle (Mie) scattering, Tp; and, selective gaseous absorption , t . Hence, the total optical 
depth is given as ® 

(12) 

The molecular- scattering optical-depth is proportional to A . This accounts for the blue 
sky, since the short wavelengths are scattered much more than the longer wavelengths. The 
wavelength-dependence of particle scattering is proportional toA^ where is dependent upon 
particle size. For extremely small particles, «x approaches 4.0 (as in ikiylelgh scattering) 
and for large particles (clouds, fog, etc. )<?< approaches 0.0, which Indicates gray scattering. 
The gaseous absorption optical depth is highly wavelength-selective, depending on the location 
of the various absorption bands of atmosirtieric constituents. Figure 28 illustrates the gener- 
al nature of atmospheric beam transmittance of solar radiation. 

Total Incoming Solar Irradiance 

The total solar irradiance, H, Is the sum of the direct-beam solar irradiance, on a 
horizontal surface, and the diffuse skylight, S; hence, 

H = I cosG. + S (13) 

where I is the normal incident direct beam. This normal incident direct beam can be given, 
in terms of atmospheric effects, as 

I = !(, (14) 

where I* is the extraterrestrial solar irradiance, commonly referred to as the "solar constant". 
For clear" days the diffuse skylight is normally blueish, which is indicative of the domi- 
nance of molecular (Rayleigh) scattering. However, the magnitude and spectral distribution 
of the diffuse skylight is a complex function of particle size, density, solar evelatlon and 
azlmui^.i, optical depth, ground albedo, particle scattering phase function, particle albedo, 
and index of refraction. For representative measurements of total and diffuse irradiance 
see reference 22. 

The diffuse skylight is also strongly influenced by clouds. Generally^ clouds increase 
the skylight and change its spectral distribution from blue to white. Clouds can also strong- 
ly reflect solar radiation and cause "bright spots" (as opposed to shadows) on the ground, 
which can very significantly affect remote sensing of ground targets (ref. 23). 

The total Incident solar irradiance is very strongly influenced by the solar zenith angle 
cos G* (equation 13 ). Since the solar zenith angle is that angle between the target's 
normal and its planar surface, topography plays a significant role. For example mountainous 
terrane has very complex topography comprising north-, south-, east-, and west-facing slopes 
of various inclinations. 

Atmospheric Rath Radiance 

Atmospheric path radiance results from molecules and large particles scattering sunlight 
in the field of view of the sensor. The magnitude and spectral distribution of path radiance 
is determined by all the variables that were given for diffuse skylight. 

Atmospheric path radiance presents a special problem because it cannot be directly meas- 
ured with ground-based instruments. Currently there are three techniques for deriving path 
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radiance : 

(1) Calculations using radiative transfer models 

(2) indirect ground-basod measurements 

(3) combination of model calculations and ground-based measurements 

Radiative- transfer models have been the subject of consideirable research. A review and dis- 
cussion of such models is given in reference 12. Models are useful for parametric studies 
of the influence of the atmosphere, given hypothetical situations. However, for specific 
conditions they are difficult to apply because, without measurements, one simply cannot con- 
fidently choose a "representative” atmospheric state. This is due to the fact that atmos- 
pheric aerosols and water vapor are known to be extremely variable as a function of geographi- 
cal locations, altitude, and local conditions. 

Indirect ground-based measurements of path radiance have been employed by Duntley (ref. 
2U) and Rogers, (ref. 25). This technique uses ground-based measurements of downward scatter- 
ed sky radiance at a geometry (sun-sensor) similar to the satellite- sun- target geometry. The 
measured sky radiance is then extrapolated to yield a path radiance for the satellite sensor 
(see ref.*s 2U euid 25). 

The third technique was used recently by Hulstrom (ref. 22) to analyze the performance 
of the Skylab EREP sensors (S-I90 A, S-I9I, and S-192). This technique uses ground-based 
measurements of optical depth/transmittance, target reflectivity, ground albedo, and total 
solar irradiance (see equation lO); then, the optical depth and ground albedo measurements 
are used as inputs for a radiative-transfer computer model that calculates path radiance. 

This technique eliminates uncertainties in atmospheric state, optical depth, and total (direct 
plus diffuse) irradiance. However, values for the particle scattering, phase function, and 
albedo have to be assumed. 

A fourth technique of deriving path radiance is being attempted in this program. This 
technique employs the satellite data itself and helicopter measurements of selected ground 
targets. If one plots equation 10, as shown in Figure 29, it can be seen that the y intercept 
is the path radiance and the slope is equal to the product HT/ir . The natural target reflect- 
ivities are measured from a helicopter platform. This is performed by having nearly identical 
spectral radiometers measure the target radiance, , from the helicopter; and, the other 
radiometer measure the total incoming solar iriradiance, H. The target reflectivity can then 
be dst irmined by 

P _tr_N^ (assumes Lambertian target) (15) 

Hence, one can determine the path radiance and the prodtjct HT/.j- . Therefore the reflectivity 
signatures of other (unknown) targets can be derived by 

P "" xtr/HT (16) 

The advantages of this technique are 

(1) it measures path radiance as "seen" from the satellite sensor 

(2) it measures the combined atmospheric effects of total solar irradiance, H, and atmos- 
pheric transmittance, T. 

(3) if the selected natural ground target reflectivities remain nearly constant over 
periods of time (as was found to be valid for this test site) numerous simultaneous 
ground-based measurements are not required in order to calibrate atmospheric effects. 

Tiie possible short-comings of this technique are 

(1) It assumes fairly unifom atnx;spherlc conditions of the surface area within which 
the natural target "call brat ioii' sites are located. 


(2) The accuracy of the technique depends on the accuracy and repeatibility of the 
measurements of natural target reflectivity, and radiances measured by sensors 
in the satellite. 

(3) It assumes that the sun angle with respect to the targets, both natural "calibratior 
and unknown targets, are similar. This is apparent when one considers the sun angle 
dependence, I cos Ooof the total solar irradiance and the product HT/fl- . 

Initially, during the summer and fall of 1975 » this technique will be employed for lANDSAT 
data. 


Pour techniques of deriving atmospheric path radiance and atmospheric effects have been 
presented. Each has it's advantages and short-comings. Unfortunately almost no data or 
experiments exist for checking and comparing the accuracy and validity of these techniques. 
All four techniques will be attempted during this program, and such comparisons subsequently 
made. 


FIELD MEASUREMENTS OF ATMOSPHERIC EFFECTS 
AND TARGET REFLECTIVITY 

Brief descriptions of the techniques and Itibrumentation used to measure atmospheric 
effects and target reflectivity are given. For a more detailed description see Hulstrom 
(ref. 22). 


Atmospheric Effects 

Atmospheric effects that are directly measureable by ground-based instruments are the 
total solar irradiemce, H, and the beeun transmittance, T. 

The total solar irradiance was measured with an I.S.C.O. (instrumentation Specialties 
Coopany) model SR (spectroradiometer) equipped with a SRR ( spertroradiometer recorder) unit. 
This instrument consists of the SR and SRR units and a 6 ft fiber-optics bundle that has a 
flat teflon diffuser, l80“ F.O.V. , aperture at the end of the bundle. Light travels through 
the diffusser, through the fiber-optics bundle into the SR unit, where it is passed through 
a chopper, slit, and monochromator. The monochromator is a wedge interference filter. The 
wavelength of the instrument is varied by lengthwise movement of the wedge interference filter 
between the slit and the sensor. The wavelength region from 38o to 750 nm is measured with 
a silicon ~junct ion i>hotocell; and, the region from 750 to 1350 nm is measured with a germanium 
Junction photocell. The 1/2 bandwidth is 15 nm from 38O to 750 nm, and 30 nm from 750 to 
1350 nm. The total solar irradiance was measured by positioning the flat teflon diffuser so 
that it was horizontal/level. Hence, the 180®F.0.V. diffuser was exposed to the total irradi- 
ance (direct + diffuse). The diffuse skylight can be measured by simply shading the diffuser 
from the direct sunlight. A complete spectral scan, from 3^0 to 1350 nm, requires about 3 
minutes. During this time interval, the total broadband irradiance (400-1100 nm) was continu- 
ously monitored with a Y.S.I. (Yellow Springs Instniraent Co.) model 68 pyranometer. This was 
done in order to assure that stable conditions existed during the I.S.C. 0. spectral scan. 

The complete I.S.C.O. unit was calibrated with an I.S.C.O. S.R.C. (Spectroradiometer 
Calibrator) unit, which uses a ribbon-filament tungsten lamp, and a N.B.S., coil filament type 
lamp (1963 standard No. E.P. I. 1523). These standards were used to derive an absolute cali- 
bration factor, C.F. , for the I.S.C.O. unit. This is derived by viewing the standard Irradi- 
ance, Iq_t.» and recording a raw intensity I . The calibration factor is then 

C.F. - ^STD/I^ (17) 

The raw field measurpment, I^, Is then multiplied by the calibration factor In order to ob- 
tain absolute quantities. However, this calibration factor Is only valid for measurements of 
normal-incident (to the (llffuser) Irradiance. For measurements of total solar irradiance, 
where the sur Is at various zenith angles (angles off of the normal to the diffuser). 
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lei source as measure-i at a Co* ancle tc tir - * Intensity of a constant paral- 

.eusurel on nonsal incidence (0*)?^InVre c:nerTr\eri“ >’ 


I(0)=I (^~o) cos 0 


( 18 ) 
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extensive measurements of collimated light were“,mdfor ' '•®=P°"S® »f the I.S.C.O.dlffuser, 
resulted In a cosine correction factor, ^C. C. F. , ^ ^ incident angles. This 


C.C.F. ^ I (0=0 ) cos 0 
I (0) 
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is made. Thf J-^Intr^cepr is lo , l^nc^vr cL'"brcS^^^^^^^ relative airmass-sec 0 , 
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measurements of direct solar irradiance the TSro'^mU instrument. Fur the 

This resulted In beam transmittance ?or^^ l3!, ' uiulpped with a 3’ collimator, 

example Is shown In fl^Sc 3l! Xre 

Within and near the Colorado tosi site ^o^^r sites located 
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atmosphere was fairly uniform. ^ shown in Figure 31, the 
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to measure target reflectivity. 
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For ground measurements of target reflectivity, standard "gray" and standard ' ’’white" 
reflectance panels are used to determine target reflectivity. These panels are accurately 
calibrated with a Beckman DK-2A/Cier-Dunkle Integrating Sphere. Field measurements consist 
of the radiance of the "gray" card, N , the radiance of the target, Nt , and the radiance 
of the ’Vhlte" card, The target ^Reflectivity ceui then be derived by 


and 



= X pgc 

gc 


(23) 


p = ^ X P»C 
^ wc 


(2U) 


where p and P are the reflectivities of the "gray" and "white" cards. The two derived 
target rerlectivlties are then averaged. This technique was verified by plotting radiance, 

N and , versus reflectivity, P and P , and extrapolating this plot to zero radiance 
ana reflectivity. If the card c^iorations 8Fe accurate and if the field measurements are 
accurate, the plot should go through the origin. Several of these verifications revealed 
the technique and measurements to be consistent to within 1 to 2^ reflectivity. Examples 
of absolute reflectivity for 15 targets and their backgrounds are listed below; 


Target Description 


Identification 

Number 


Pikes Peak Granite 
Lichen Rock Cover 
Soil/Grass Background for 
Pikes Peak Granite 
Cripple Creek Granite 
Soil/Grass Background for 
Cripple Creek Granite 
Volcanic Fine-Grain Mafic 
Soil/Grass Background for 
Volcanic Mafic 
Meadow Grass 
Volcanic Crystal Tuff 
Yellow Grass Background 
Volcanic Andesite 
Background Soil for Vol- 
canic Andesite 
Background Grass for Vol- 
canic Andesite 
Background Soil at Cameron 
Mountain 

Background Soil/Grass at 
Cameron Mountain 


1 

2 

3 

U 


5 

6 

7 

8 
9 

10 

11 
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13 

lU 

15 


ORIGINAL PAGE lb 

rgjt or POOR QUAUTVl 

No. Percent Reflectance 



0.5 to 

0.6 to 

0.7 to 

0.8 to 


0.6 

0.7 

0.8 

1.1 1 

1 

lU.o 

18.U 

22.6 

25.6 

2 

23.6 

28.5 

36.1 

35.2 

3 

12.7 

lU.8 

20.7 

23.4 

U 

17.3 

22.6 

24.5 

25.6 

5 

11.1 

13.6 

19.3 

22.6 

6 

11.7 

14.5 

19.1 

19.4 

7 

10.8 

16.8 

19.5 

26.0 

8 

5.8 

5.4 

30.8 

42.5 

9 

27.5 

31.8 

35.8 

4o.4 

10 

11.8 

13.6 

19.1 

18.9 

11 

lU.5 

17 7 

19.8 

20.9 

12 

9.3 

13.: 

18.0 

18.7 

13 

6.5 

8.5 

12.7 

17.3 

l4 

10.8 

i 6.3 

19.2 

23.2 

15 

9.6 

12.8 

19.8 

24.1 


In order to measure signatures from mixtures and for areas that are on the same order 
of size as the satellite sensors, helicopter measurements are employed. Helicopter measure- 
ments consist of using one EXOTECH radiometer, with a 1“ FOV, in the helicopter to measure 
the target radiance, and one EXOTECH radiometer located on the ground to measure total 
incoming solar irradiance, H. The reflectivity, (> is then 

p (25) 

In addition, the helicopter measurements of will be compared with the computer simulations 
of the radiance signature, 

SUMMARY AND CONCLUSIONS 

Methods have been discussed for simulating spectral signatures of mixtures of two terrain 
classes for automatic analysis of LANDSAT-1 MSS data using on-site measurements, simulated 
LANDSAT data, and actual LANDSAT-1 data. Tests of the methods for estimating component-class 
spectrcd. signatures and simulating two component- class mixture signatures with simulated 
MSS data indicated that improvements in classification results over conventional con^xjnent 
class analysis are possible using those techniques. Applying these techniques to actual LALD- 
SAT-1 I^S data of wildland areas showed an increase in classification information over conven- 
tional analysis, but no quantitative accuracy analysis could be made due to difficulties of 
estimating mixture proportions from the control or "ground truth data. 

It was found that mi solassifi cation of pixels as mixtures can occur when the simulated 
spectral signatures appreximato those cf component classes or other mixtures in the scene. 
There was also evidence of misclassificatiors due to mere than two component classes existing 
within a single pixel. Within-clnss variability, slope and aspect variability, and the sen- 
sitivity of the scanner in detecting changes in mixture proportions may also contribute to 
degrading the classification porformance using these methods. 


These techniques may pi*ov"‘ valuable In areas where vegetation masks the characteristic 
spectral response of the underlying geologic material , and future work should include an ex- 
amination of this application, amphasis should also be placed on methods of acquiring better 
"ground truth"information upon which the spectral signature estimation and simulation techni- 
ques are based. 


I^pertles of the atmosphere play an Important role in affecting the spectral 'radiance 
received at the satellite by the scanner. Some atmospheric properties can be measured on 
site. Others must be determined by comparing satellite data vlth the on-slte measurements. 

Once their signatures have been determined, it may be possible for some large natural terrain 
features to serve as calibration panels from which atmospheric properties can be calculated 
by the satellite data. When the desired remote sensor product is an image having only rough 
spectral Information, atmospheric effects caui be ignored. However, when quantitative spectral 
infonnatlon is required, atmospheric effects cannot be neglected. 
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Figure 1. -location of test site. EMC=Eleveninile Oanyon Reservoir Study 
Area, P= Pikes Peak, C^Cripple Creek, COOanon City. 



EXTRACTED SPECTRAL SIGNATURE PLOTS 



Flfiur« 2, -Spectral responie of two component terrain classes (A.B) 
and a mixture of them (C) extracted from LANDSAT-1 MSS data. 
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Figure 3* "Comparison of total solar Irradiance and reflectivity of coniferous 
forest and Pikes Peak Granite. 
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Figure 5.-Spectral radiance of Pikes Peak Granite, coniferous forest, 

tures of the two. Curves, date, time, and other documentation 
in fig, U. 
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Figures 6 and 7. -Spectral inherent contrast of Pikes Pt*ak Granite (C) target 
with background of Coniferous Forest (F), and mixtures of granite and 
forest. Curves: G F, G F, ?»7(V)t.O ^ F, U=so% q 

F, 5-KKVit F. 
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EXTRACTED SPECTRAL SIGNATURES VS. CALCULATED MIXTURE SIGNATURE 



^ 567 


MSS BAND 


Figure 8, -Spectral response of two component terrain classes (A,B) 
and a mixture of them (C) extracted from LANDSAT-1 MSS data, 
and the mixture response (D) calculated using equation U. 

A, B, and C, same data as in fig. 2. 




looi r 

looi G 

looit w 

75lt G 25^ F 
50lt G 5016 P 
25')t G 75^ F 
75^ G 25i W 
50 ^ G 50^t W 
?5^t G '^•5^ W 


Figure 9 . -Microfilm graymap of 
band represents a field 
labelled. F=Coniferous 


9-class simulated data set. Each horizontal 
of simulated data for a given class, as 
forest, G=grassland, VMrater. 


u 
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Figure 10. -Mean spectral response for two compone^ ^sses and Jjod^ed mix- 
ture classes. Curves: G-100% grassland, F-100^ devStlins 

G 50 % F, E- 25 ^ G 75% F. Dotted lines represent standard deviation . 





100 ^ 0 
G25i r 

5016 G 50 % F 
25 % G 75% F 
100 % F 


FlKure ll.-Microfilm classification display for component-class 
analysis. Blanks denote unclassified points. Symbol 
are:‘ • =grassland(G) , .=forest(F). C,D,E same as In fig. 1 . 



£sns“.“r.‘"r“" 

classes and mixtures as in fig. 10. 




RESPORSE 



Figure 13, -Mean spectral response curves for simulated Omssland- 
Water mixtures. Dotted lines indicate standard deviations. 
Letter symbols are same as in matrix on p. 12. 



SIMULATED SPECTRAL SIGNATURE PLOTS 



MSS BAND 


Figure lU.'Mean spectral response curves for component and 
mixture class signatures used for 9*class analysis. 
Dotted lines Indicate standard deviations. Letter 
symbols are same as In matrix on p. 12. 
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Figure 15. -Microfilm classification display of 9-class analysis. True class 
field symbols noted at left of display are same as in matrix on p. 12. 



Figure l6. -Microfilm classification display for 9-class analysis with calculated com- 
ponent class signatures from uniform data. True class field symbols noted at 
left of display are same as in matrix on p. 12. 
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Fl«ir« l8.-Mlcrofll« ciaa«ification displays for conponent class a^ysls 

mated mt^ n vectors and covariance matrices, a) Spectral signatures estimated 

trtm onlfonn trilnlog set., b) Spectral ^„rfirtd 

training sets. Symbols are: ■ -Grassland, \ -Forest, ^ ^ 

Classes noted at left of displays have same symbols *s In matrix on p. 
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Figure 19**Microfllin classification display of 9 ~c1ass analysis using 
calculated component class mean vectors and common covariance 
matrix. True class field symbols noted to left of display are 
same as in matrix on p. 12. 
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Figure 21.-Mlcroflljn claeslficutlon display for 9-da«« analysis using cal- 
culated mean vectors and averaged cosnon covariance matrix. True class 
fields noted at left of display are same as in matrix on p. 12. 


Figure 20.-Hicrofilm classification display of 9-class analysis using cal- 
culated mean vectors and cominon covariance matrix. True field class 
symbols noted at left of display are same as in matrix on p. 12. 





Figure 22.-MicrofilJtt graynap of MSS Band 5 of the Elevenmlle 
Canyon Reservoir Study Area. 
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-Mlcrofil* classification display 

Werwlr Te.t Are. vUK 5 

.-Fbrest, V -Grassland, 7 ^ 

, .Wet Meadoir, Blank Unclassified. 


of Elevenalle Canyon 
Classes. Synbols are: 
-Mountain Grassland, 





classification map of Qevenmlle Canyon Reservoir study area with 
classes and 9 mixture classes. Because of limitations of the mlcro- 
fila display, separate symbols could not be shown for the mixtures- however for 

”” ‘'y » <»n».ntion.l llne-prlrter. * 


•Grassland (g) 

> Forest (F) 

'^Mountain Grassland (mg) 
'Wet Meadow (WM) 

• Whter (W) 


Mixture classes: 


X 


7516 G 2516 P, 501fc G. 
7516 G 2516 WM, 5016 G 
7516 F 25% MG, 50% P 


50% P, 25% G 75% F, 
50% WM, 25% G 75% WM, 
50% MO, 25% F 75% MG 


Blank, unclassified 






Figure 25. -Microfilm classification map of Elevenmile Canyon Reservoir study area with 

5 component classes and 10 mixture classes. Because of limitations of the micro- 
film display, separate symbols could not be shown for the mixtures; however, for 
study purposes they were printed separately by a conventional line printer. 
Component classes: 

• Grassland (G) 

\ Purest (r) 

\ Mountain Grassland (MG) 
vWet Meadow (WM) 

V Water (w) 

Mixture classes: 

■ Silo G F 
• Silo G 3?i W 

X - 671t WM 331t w, 331t WM 671. W, 331t G SH F, 

331t G 671 W, 671 F 331 MG, 33l F 6t1 MG 
671 G 6t1 WM, 331 G 671 WM 


Blank, unclassified 
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Figure 26. -Mean Spectral Response Curves for Msnltou Study Plots and Estioated 
Cksnponent Qasses. Symbols are: A« Estimated Forest Curve » ^Flot 

OPlot #3, I>-Plot #1, E»Plot |U, F.piot #6, G-Plot #5, H-Estimted Back' 
ground Curve. 




S(XAR IRRADIANCE 
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Figure 27. -Satellite remote sensing of target on the ground. 



WAVELENGTH, MICROMETERS 

Figure 28. -Typical spectrum of solar energy reaching the surface of the earth. 

From Handbook of Geophysics and Space Environment, S. L. Vklley, Editor 
(New York: McGraw-Hill Book Company, Inc.), p. I0.2 


RADIANCE RECORDED AT SATOXITE 
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COSIHE CORRECTICW FACTOR 


2.0 




Figure 30, -Cosine correction Actors for TSno 

near-infrared (b) spec^JL t f r 0 “d 

optics probe, in b each wiw* * * * using six-foot fiber 

N.M. (top) to 750 N.*M. (bottom) /*^*^***"^“ * wavelength from 1350 
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Villa Grove 

-■Transmittance for Relative Air Msss, 
m* 1.0 

a 

Hillside 


Granite Hills 


© 

Hartsel 


WAVELENGTH, NM 


Figure ll^'Fractlonal transmittance for sites in south-central Colorado. All sites are separated by about 
80 km. Granite Hills site lies within study area of figure 1. All data from October 2, 197U. 






ORIGINAE PAGBH 
OF POOR QUAUni 



FIGURE 32. Color com|)osUe LANDSAT-1 image of test site region. Image 
1028-17135 of August 20, 1972. Annotations same as for figure 1: 
P=Pikes Peak, C=Cripple Creek, CC=Canon City, EMC=Elevenmi le Canyon 
Reservoir. Because rectangular pixels are printed as squares by 
the computer the image is distorted, both in azimuths and distance. 
The scale shown is an average. 
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FIGURE 33 (ON FOLLOWING PAGE) . -Computer-generated recognition map of 
test site. Compare with figure 32 for location of features. Made 
by Environmental Research Institute of Michigan (ERIM) from computer 
compatible tapes from LANDSAT-1 image 1028-17135 of August 20, 1972. 
Because rectangular pixels are printed as squares by the computer, 
the image is distorted both in azimuths and distances. The scale 
shown is an average. 

EXPLANATION OF COLOR CODE 


1. Niobrara Shale 

2. Pierre Shale 

3. Limestone 

4. Canopy of Ponderosa Pine greater than 25 percent 

5. Dakota Sandstone 

6. Fountain Formation 

7. Volcanic and plutonic rocks, undivided 

8. Granite-Tundra composite 

9. Granitic grus 

10. Meadow 

11. Water 

12. Cloud shadow 

13. Cloud 




APPENDIX A 


RECOG Computer Program Blocks 

At Colorado State University, conventional automatic analysis of MSS data is accomplished 
through the use of a series of pattern recognition programs called RBCOG. RECOG consists of 
six program blocks %rhich are a modification of an original version, called LARSYS, devtioped 
at Rxrdue University (ref. 10). These programs provide a logical procedure for processing 
MSS data using supervised learning techniques. A brief description of each program block 
or phase is presented here. 

Riase 1 


Ihase 1 is a display routine that provides the user with a computer line-printer or 
microfilm representation of the scene that is to be aut<anatically analyzed. The MSS data are 
disi>layed as a grayinap representing the radiation response from each pixel. A range of radi- 
ation responses (either specified or default) are coded as a symbol and displayed for a select- 
ed wavelength band. This display provides a pictorial representation of the MSS data from 
which the user may delineate field boundaries of terrain classes that he may wish to map. 

Riase 2 


The fields identified on the Phase 1 graymaps are used as training sets from which the 
mean spectral response and standard deviation vectors and correlation and covariance matrices 
are determined by Phase 2. The mean spectred response vector and covariance matrix provide 
a statistical spectral signature for a terrain class that is used in a later Phase to auto- 
matically classify each point in the MSS data set. 

Phase 3 


When the multispectral scanner has a large number of channels available, processing the 
data using all of the information becomes quite expensive. Phase 3 is designed to select 
a subset of optimum channels for identifying all of the terrain classes utilizing divergence 
criteria. 

Phase 4 


The spectral signatures obtained for each terrain class from Phase 2 can be analyzed as 
to how well they represent the class by selecting a subset of the MSS data and classifying 
it with Ihase U. Ihase h is designed as an instructional mode and allows classification of 
the data set with three algorithms: LEVELS, a level-slicing routine; EUCLID, a Euclidian- 

distance routine; and GLIKE, a maximum-likelihood routine. GLIKE is the algorithm used to 
classify the data in the next Phase so it is valid to test the representativeness of the 
signature set with it. This allows refinement of each spectral signature by redefining the 
training set to discrad any point that would tend to make a signature unrepresentative of its 
respective terrain class. 

Phase 5 


Phase 5 is the actual classification mode for the RECOO pattern recognition sequence. 

The mean vector and covariance matrix for each terrain class are used with the likeli- 

hood decision rule GLIKE to classify each pixel in the specified MSS data set. As the classi- 
fication results are generated they are written onto a magnetic tape and stored as a permanent 
file. 

Phase 6 

The final step in the processing scheme. Phase 6, displays the results generated by Phase 
5. The user is given the option to specify a threshold level which sets a confidence limit 
for the classified data points. This is designed to eliminate false classifications of data 


which do not fit any of the terrain classes. Current display modes available with l^ase 6 
include a thematic map produced on the computer line-printer and/or microfilm similar to 
those produced by Riase 1, except each pixel is identified by a symbol or intensity level 
(line printer and microfilm respectively) representing a terrain class. A more detailed de- 
scription of RECOG can be found in ref. k and 13. 

APPENDIX B 


EQUATIONS 

The meem vector and covariance matrix that describe the spectral response of a terrain 
class are determined from a set of MSS data points known to contain that class. The statisti- 
cal equations used to compute them are discussed below. 


Since the radiation response from a ground resolution element sensed by a multispectral 
scanner can be described by the column vector 



(26) 


where each x comi»nent ^presents the radiance recorded in a spectral channel, we can find 
expressions for m' and C: in terms of the radiation responses X^. (ref. 10) 

The mean vector is given by the column vector 



(27) 


where m^^ is the mean spectrcd response in wavelength band i (for ERTS-1 i=l ,U) given by 

n^ 

Xik (28) 

k-1 

where n^ is the number of pixels (sample points) in the training set describing terrain class 
A, and k is the sasple point index. 


The covariance matrix which indicates hew the radiation response in one MSS channel varies 
with the response in the other channels can be described as 


• • • 


where f. Is the covariance between channels i and j given by 


Ai ■ "j’ * 

k=l 

Bie steuidard deviation for channel i in class A is , and the correlation coefficient 
between channel i and channel j for class A is given by 
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